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Abstract. The hypothesis that the Damped Lya systems (DLAs) are large, 
galactic disks (Milky Way sized) is tested by confronting predictions of mod- 
els of the formation and evolution of (large) disk galaxies with observations, 
in particular the Zinc abundance distribution with neutral hydrogen col- 
umn density found for DLAs. A pronounced mismatch is found strongly 
hinting that the majority of DLAs may not be large, galactic disks. 



1. Introduction 

It has been proposed that the Damped Lya (DLA) systems seen as broad 
absorption lines with neutral hydrogen column density N{HI) ^ 2 10^" 
cm~-^ in quasar spectra are large and massive (Milky Way sized), galactic 
disks with maximum circular speed Vc <^ 200 km/s (e.g., Prochaska & Wolfe 
1997, 1998). If true, this is of profound interest in relation to theories of 
the formation and evolution of disk galaxies, but it is obviously important 
to test this hypothesis in as many ways as possible. In this contribution I 
present a test alternative to the ones used by Prochasca & Wolfe. My test 
is based on fairly elaborate and realistic models of the formation and evo- 
lution (including chemical evolution) of large and massive disk galaxies, in 
particular on comparing the observed distribution of Zinc abundance with 
neutral hydrogen column density found for DLAs with model predictions. 

In section 2 I briefly describe the models, in section 3 the data, in section 
4 the model predictions are confronted with the data and finally section 5 
constitutes a brief conclusion. 
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2. The models 

The star-formation rate is assumed to depend on gas surface density as 



OC w{Egas) , (1) 



dt 

where N=l-2, and the weight factor 

wCSgas) = 1 (2) 

for models without a star-formation threshold and 

~ , < Et// 

wi^gas) = I (3) 

for models with a star- formation threshold given by (e.g., Kennicutt 1989) 

^TH = a ^ 11 (-^) (— ) (-^) Mo/pc2 , (4) 

3.36 G 0.7 Kq 6km/s 

where k is the epicyclic frequency 

c is the velocity dispersion of the gas, a is a dimensionless constant of the 
order unity and k© is the epicyclic frequency at the solar distance from the 
center of the Galaxy. I have adopted a = 0.7 and c = 6 km/s (Kennicutt 
1989). Infall of primordial gas, but no gas out-flows arc included in these 
models of large and massive disk galaxies. Chemical evolution is calculated 
assuming a Scalo (1986) initial mass function (IMF). Oxygen is assumed to 
be produced in massive stars only (M > 8 M©) and instantaneously recy- 
cled by Type II supernova (SN) explosions. Iron is assumed to be produced 
partly by massive stars {M > 8 M©) and instantaneously recycled by Type 
II SN explosions and partly by less massive stars and recycled by Type la 
SN explosions with a timc-dclay of tsNia = 1 Gyr. Unprocessed gas, locked 
up in stars, is assumed to be non-instantaneously recycled. 

The proportionality constants in the star-formation laws (eq. [1]) are 
determined by applying eq. [1] to the local, solar cylinder and fitting data 
on the local gas fraction, G-dwarf distribution etc. 

The global properties of the large disks as a function of time (or red- 
shift) are determined by assuming that an average, spin-parameter A = 
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Detections, all z, <z> = 1 .91 
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Figure 1. The observed distribution of [Zn/ H]dla vs. neutral hydrogen column density 
N{HI) for the 18 detections in the combined Pettini et al. sample. 



0.05, galactic disk of either early or late type (see below) has an exponen- 
tially decreasing stellar profile with radius R in the disk (e.g.. Freeman 
1970) truncated at 4 stellar scale-lengths (e.g., van dcr Kruit 1987) and an 
exponentially decreasing oxygen gas abundance profile with radius in the 
disk (e.g., Zaritsky et al. 1994, Garnett et al. 1997). Early type disk galaxies 
(Sab/Sbc) arc assumed to have an average star-formation history parame- 
ter b, defined as the ratio of the current to the past average star-formation 
rate, of 6 ~ 0.33 and late type disks (Scd/Sdm) to have an average b ~ 
1.0 (Kennicutt et al. 1994). An Ellis (1983) present day iz=0) morpholog- 
ical mix of E/SO:Sab/Sbc:Scd/Sdm = 0.28:0.47:0.25 is assumed and the 
models are averaged over spin-parameter A. The models thereby cover the 
full range of morphology and global surface brightness of large, galactic 
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Upper limits, all z, <z> = 2.56 
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Figure 2. Same as in Figure 1, but for the 16 upper limits 



disks and are fully constrained: There are no free parameters left, except 
for the star-formation law index N and the absence or presence of a star- 
formation threshold and the conclusions reached in this work turn out to 
be quite insensitive to this for reasonable values of N {N=l-2) - see also 
below. 

3. The data 

Observational data on 34 DLA systems (including Zinc abundances or lim- 
its) have been taken from Pettini et al. (1994, 1997). Zinc is used because 
it is not or only very slightly affected by dust depletion. Figure 1 shows 
[Zn/H] versus neutral hydrogen gas column densities N{HI) for the 18 
DLAs with Zinc detected. They have a mean redshift of < z >=1.91. In- 



ARE DAMPED LYa SYSTEMS LARGE, GALACTIC DISKS ? 



5 



Stellar mass return Included using a Scale 1986 IMF 
n=0.3, z=1.91, N=1.5, dust-corr., Mix, X & proj. area weighted 
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Figure 3. Predicted DLA abundance distributions at different N{HI) assuming 
that DLAs arc large, galactic disks similar to the Milky Way: For models with no 
star-formation threshold. Dotted line histogram: N{HI) = 2 — 10 10^° cm~^, Solid line 
histogram: N{HI) = 10-50 10^" cm'^ 



terestingly, higher Zinc abundance corresponds to lower N(HI), which is 
not to be expected if the DLAs are large, galactic disks, even when effects 
of dust obscuration of quasars are taken into account - see below. Figure 
2 shows [Zn/H] versus N{HI) for the 16 DLAs with upper limits on the 
Zinc abundance. The mean redshift of these is < z >=2.56. 



4. Models versus data 



Figure 3 shows the model predicted [Fe/H] distributions for N{HI) = 2- 
10 10^0 cm-2 and N{HI) = 10-50 10^° cm-^ for the no star-formation 
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Stellar mass return Included using a Scalo 1986 IMF 
0=0.3, z=1.91, N=1.5(TH), dust-corr.. Mix, X & proj. area weighted 
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Figure 4- Same as in Figure 3, but for models with a star-formation threshold 



threshold models with star-formation law index A'"=1.5. The distributions 

have been weighted over morphological mix, spin-parameter (and hence 
effectively global surface brightness) and projected area and corrected for 
the effects of dust obscuration of quasars using the clever methodology 
developed by Fall Sz Pei (1993). Figure 4 shows the same thing, but for the 
models with a star-formation threshold (and still iV=1.5). 

Neither of the two sets of model distributions appear to match the data 
at all (assuming that [Zn/H] ~ [Fe/H] as is the case for Galactic stars 
of all abundances). More quantitatively, Komolgorov-Smirnov (KS) tests 
show that both sets of model distributions can be ruled out as fits to the 
data with more than 99% confidence (this result is more generally obtained 
for all values N=l-2 and also for Wyse &; Silk (1989) type star-formation 
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dN/dz for DLAs with Zinc detection 

Data: APM survey based; Models: Disks with V^S200 km/s, N=1.5 
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Figure 5. Predicted number of DLAs per rcdshift interval for the models with no 
star-formation threshold - see text for details. Points with errorbars: Data, Lines: Models. 



laws). 

The mismatch between the data and the models is further illustrated by 
Figure 5. The figure shows the number of predicted DLA systems per red- 
shift interval with Zinc detected, dN/dz, for N{HI) = 2-10 10^° cm'^ and 
N{HI) > 10^^ cm~^ versus redshift z. The models have been normalized 
using observational information about the luminosity function and sizes of 
present day (z=0), large and massive disk galaxies (with Vc <^ 200 km/s). 
The model predictions have been weighted and corrected as described above 
and furthermore corrected for the effect of a Zinc detection threshold as a 
function of N{HI). Figure 5 is for the models with star- formation law index 
A^=1.5 and no star- formation threshold. Also shown are APM survey based 
observational results of Storrie-Lombardi et al. (1996) corrected by an ob- 
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servationally deduced Zinc detection probability as a function of redshift 
z and N{HI) derived from the Pettini et al. (1994, 1997) data. As can be 
seen from the figures the two sets of model predictions lie more than two 
orders of magnitude below the observations - the situation is even worse 
for the models with a star-formation threshold (not shown). 

5. Conclusion 

The analysis presented above suggests that the DLA systems may not in 
general be accounted for by large and massive disk galaxies with Vc ^ 200 
km/s, not even when account is taken for the existence of some very large 
and low surface brightness disk galaxies (physically: high A galactic disks). 
The true situation is most likely quite complicated, but it is tempting to 
associate the DLAs at low redshift with the large population of faint, blue 
dwarfs and at higher redshift with the large number of smaller galaxies 
expected at such redshifts {z ~ 2-4) in hierarchical galaxy formation sce- 
narios (e.g., Kauffmann 1996, Mo et al. 1998). DLAs with indications of 
large internal velocities (~ 200 km/s) may, at least partly, be dynamically 
unrelaxed systems (e.g., Haehnelt et al. 1998). 
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